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^SJ ■ Abstract 

A fluid, like a quark-gluon plasma, may possess degrees of freedom indexed by a group 
variable, which retains its identity even in the fluid/continuum description. Conventional 
qq ■ Eulerian fluid mechanics is extended to encompass this possibility. 

<n : 
o ■ 

^ '. 1 Introduction 

o : 

The dynamics of fluids is described by a classical field theory, whose origins lie in the nine- 
teenth century, like Maxwell electrodynamics with which it shares some antecedents. The 
electromagnetic theory has enjoyed much development: it passed into quantum physics, and 
the resulting quantum electrodynamics served as a model for its non-Abelian generalization, 
Yang-Mills theory, which today is at the center of a quantum field theoretic description for 
fundamental physics. We believe that fluid dynamics can undergo a similar evolution and 
^ ■ thereby be generalized to fluid systems with non-Abelian charges. 

The quintessential example of a physical system with non-Abelian charges is the quark- 
gluon plasma. High energy collisions of heavy nuclei can produce a plasma state of quarks and 
gluons. This new state of matter has recently been of great interest both theoretically and in 
experiments at the RHIC facility and at CERN. In fact, there is growing evidence that such a 
state has already been achieved at the RHIC facility [Q. 

In attempting a theoretical description, there are basically two approaches that we can 
use. Since the plasma is at high temperatures, one can argue that the average energy per 
particle is high enough to justify the use of perturbative Quantum Chromodynamics by virtue 
of asymptotic freedom. However, it is known that because of infrared divergences various 
resummations, such as summing hard thermal loop contributions, have to be done before a 
perturbative expansion with control of the infrared degrees of freedom can be set up One 
has to address also the question of chromomagnetic screening, because unlike the Abelian 
plasma, there can be spatial screening of magnetic type interactions j3|- The expected end 
result is then a good description valid at high temperatures and for plasma states that are 
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not too far from equilibrium, since one is perturbing around the equilibrium state. An alter- 
native approach, which may be more suitable for nondilute plasmas or for situations far from 
equilibrium, would be to use a fluid mechanical description. 

We begin by recalling the conventional formulation of Eulerian fluid mechanics (in the 
Abelian case), for perfect, isentropic fluids [I] The variables are the matter density field p(t,r) 
and the velocity field v(t, r); the two are united in the current density j = vp. 

The equations satisfied by these in the non-relativistic situation are the continuity equation 

d d 

- /O +V-(vp) = = -p + V-j, (1.1) 

and the Euler equation, 

d 

— v + v • Vv = force density. (1-2) 

The first expresses fluid matter conservation; the second describes forces. The force density 
on the right could be derived from a pressure function, P(p), 

force density = — VP(p), (1.3a) 

or more generally, from an enthalpy function V(p), 

force density = — W'(p), (1.3b) 

or in case of magneto hydrodynamics from the Lorentz force law in the presence of electric E 
and magnetic B fields 

force density = e(E + v x B). (1.3c) 

(Here e is the charge; the velocity of light is scaled to unity; a dash, as in (|1.3b|) . denotes 
derivative with respect to argument.) In the first two cases it is assumed that the dependence 
of P or V on p is known (= equation of state) and the system Ijl.ljl . (|1.2|) is self-contained. 
In the magneto hydrodynamical case, the the electromagnetic fields are taken to be either 
prescribed functions of (t, r), or determined self-consistently from the Maxwell equations with 
sources (p,j). 

A simplification arises when the vorticity u = V x v vanishes. Then the velocity can be 
written as the gradient of a velocity potential 9 

v = V(9 (1.4) 

and equations ()1.2j) . ()1.3b|) are integrated once to yield the Bernoulli equation 

wt e +r> = -v {P) (1 .5) 

In the relativistic case, the matter density and current density are united in a current 4- 
vector j 1 * = (p,j) = (p, pv) which is decomposed into proper density n and a proper velocity 
u^, normalized to unity. 

f = nu^ u^u^ = 1 (1.6) 

This decomposition, which entails no restriction - it is valid by definition - is called the 
Eckart factorization 5 . [Note that the Eckart factorization also holds non-relativistically 
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(pJ) = /°(l) v )) but the velocity is not normalized.] The current is conserved, satisfying a 
continuity equation, 

V = 0, (1-7) 

and the Euler force equation is presented with the help of the Eckart factorization as an 
equation for the proper velocity u. 

u^d^ u v = force terms (1.8) 

Equivalently, the non-relativistic/relativistic Euler equations can be presented as continuity 
equations for a non-relativistic/relativistic energy- momentum tensor 9^ v . 

d^ u = (1.9) 

Additionally there exists a phase space formulation, with known brackets for p and v, and a 
Hamiltonian generates the equations of motion by bracketing . Also there is a configuration 
space formulation in terms of a Lagrangian. This requires writing the velocity in terms of the 
Clebsch parameterization 6 . 

v = Vf9 + aV/3 (1.10) 

Note that in three dimensions the (p, v) phase space is 4-dimensional, hence even. Also the 
velocity vector field v involves three functions, which are captured by the three scalar Monge 
potentials (6, a, (3). [That a Lagrangian/configuration-space formulation requires presenting 
the velocity vector field in terms of potentials, which are not seen in the equations of motion 
is analogous to the Lagrangian/configuration-space formulation for a charged particle moving 
in a magnetic field: the Lagrangian involves a vector potential, while the equation of motion 
sees only its curl - the magnetic field.] 

The research problem that my collaborators and I set for ourselves concerns the question 
of how to generalize the above, when there is a variety of fluids, with each species labeled by 
a Lie group index j^, [7]. Certainly we expect that the current will inherit the group index 

— > Jo , and will be group-covariantly conserved. 

dpJZ + fabc Apt = (D^) a (1.11) 

Here A^ is a non-Abelian gauge potential, describing an external or dynamical gauge field in 
which the fluid moves. In the latter case we expect that the gauge field obeys its own equation 
of motion. 

+ fabc A\ Fr = (£>„ F^)a = J u a (1.12) 

The f a b c are structure constants for the appropriate group. 

But the question that still remains is how should the density and the velocity of the fluid 
be related to the group covariant current; viz. what are the group transformation properties 
of the density and velocity; what is the analog of the Eckart factorization ()1.6|) : which is the 
non-Abelian Euler equation that supplements the non-Abelian current conservation equation 



2 Non-Abelian Euler variables 

A possible form for the non-Abelian Euler fluid variables may be inferred from the single- 
particle equations of motion — the Wong equations [H] (over-dot denotes time differentiation) , 
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mX n _ pA4) (2.1a 

yd 



mP^ = gQ a F a P\ (2.1b) 



mQ a = -gfabcP^Qc, (2.1c) 

with A't and F% v describing a (background or dynamical) gauge field interacting with strength 
g. [We have assumed that the only active forces are Lorentz forces due to the gauge fields. 
Of course other forces can be included on the right of l)2.1bjl .] The equation of motion (|2.1cj) 
satisfied by the charge Q a is postulated to hold so that the single particle current 

J£(t, v) = fdr Q a {r) X»<j) 5{X\t) - t) <5(X(r) - r) (2.2) 

is covariantly conserved, simply as a consequence of its definition (|2.2[) . and of the equations 
(12~H) satisfied by X»,X» oc P», and Q a . 

A convenient choice of parameterization is X (r) = r, so that Ja = (p a , Jo) is given by 

Pa(*,r) = Q«(t) 5(X(*) - r), (2.3) 
J«(t,r) = Qa(t) X(t) 5(X(t) - r). (2.4) 

To arrive at the corresponding fluid equations, we first consider an assembly of N particles; 
i.e. we decorate all the co-moving variables with a particle-counting index n = (1,...,N) : 
Xn,Pn, Qan, and replace the single-particle formulas (|2.3|) and ()2.4|) with sums over n. Then 
we pass to a continuum/fluid description by letting the discrete counting variable n become 
the continuous label x, and we replace sums over n by integrals over x. 

Pa(t,r) = |A Q tt (t,x) p(t,x)-r), (2.5) 
Ja(*,r) = J (Fx Q a (i,x) X(t,x) $(X(t,x) - r) (2.6) 

The continuum version of (|2.1|) becomes (we suppress m and 5) 

X„(t,x) = Q a (t,x.)F^(t,X(t,y:))X u (t,x.), (2.7a) 



Q a (t,x) = -fabc 



A b Q (t, X(t, x)) + X(t, x) • A fe (t, X(t, x)) 



Q c (t,x). (2.7b) 



As a consequence of these, Ja of (|2.5|) . (|2.6|) is covariantly conserved. 

Let us now observe that Q2.5|) and Q2.6JI describe an Eckart factorization for the non-Abelian 
case in the following manner. In the course of the x- integration, the <5-function evaluates x at 
function %(t,r), such that 

X(t, X (t,r))=r. (2.8) 

Also there is a Jacobian factor. Thus 

p a (t,r) = Q a (t,r)p{t,r), (2.9) 
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where 



Thus we may write 



Qa(t,r) 

P(t,r) 
Ja(t,r) 
v(i,r) 



Qa(t,x) 

1/det 



x=x(*.r), 



8x3 x= X (t,r.. 

Q a (t,r) v(i,r) p(t,r) 



X(t,x) 



x=x(*,r) 



Q a (t,r) i"(t,r), 
p(*,r), v(t,r)p(t,r) 



(2.10) 

(2.11) 
(2.12) 
(2.13) 



(2.14) 
(2.15) 



Since pand j = vp are defined by integrals (|2.5j) and (|2.6|) without the charge Q a (t,x) factor, 
the Abelian current is ordinarily conserved. 



d»f = 

Together with the covariant conservation of Ja, this implies that 

f{D»Q) a = o, 



(2.16) 



(2.17) 



which we call the fluid Wong equation. Note also that the Abelian current can be factorized 
d la Eckart as in (|l.fi[) . so that the complete non- Abelian current reads 



1. 



(2.18) 



[Unlike the Eckart decomposition for the Abelian current, (|1.6|) . which entails no restriction, 
the non- Abelian decomposition (|2.18j) is specific to our assumption that a particle substratum 
underlies the fluid variables. An alternate picture is presented in Section 4.] 

These results indicate that the color degree of freedom is carried by a color density Q a , but 
all colored components move with the same velocity v or ~ a single flow characterizes the 
fluid. 

It remains to posit an Euler force equation, like (|1,8|) . which captures the dynamics. This 
we do by specifying a Lagrangian action, which will also lead to (|2.16|) and (|2.17() . Moreover 
our Lagrangian will offer the possibility of generalizing the above single flow description to 
allow for different flows in different components of the fluid 0] , [7] . 



3 Constructing the action 

The equations of motion for the non- Abelian fluid in the Euler formulation include the kine- 
matical equations: continuity Ijl.llJ) and Wong (|2,17|) that are satisfied by the non- Abelian 
current, which is factorized as in (|2.18|) . Still needed is the non- Abelian Euler force equation, 
which extends (|1,8|) . or non-relativistically, (jl.2|) . thereby specifying the dynamics. We present 
this by constructing an action whose vanishing variations reproduce the kinematical equations 
and give a model for the dynamical equation. 
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The algebra underlying the non-Abelian theory is realized with anti-Hermitian generators 
T a satisfying 

[T a , T b ] = f abc T c , (3.1) 

and normalized by 

tr T a T b = __ § ab ^3_ 2 ) 

In a canonical particle theory we expect that the algebra (|3.1j) is reproduced by Poisson brackets 
for corresponding symmetry generators. In a field theory, we expect to find a copy of (|3.1I) at 
each point in space, leading to the Poisson brackets 

{p a (r), Pb (r')} = f abc Pc (r) 6(r - r'). (3.3) 

[A common time argument is suppressed.] Upon quantization the brackets become commuta- 
tors and acquire the factor i/K. (Schwinger terms do not spoil the quantum algebra unless there 
are anomalies in the gauge symmetry. Of course, we assume that the theory is anomaly- free.) 

M 

The action which leads to the bracketing rules (|3.3[) is the field theoretic version of the 
Kirillov-Kostant form, which for a particle (not a field) reads ^U] 

where n is a normalization constant and we have taken the group SU (2) as a concrete example: 
g £ SU(2),g = exp(T a (p a ),T a = a a /2i and a a are Pauli matrices. Also one uses ()3.4(l as part 
of an action for the Wong equation of motion for a particle. 

More generally, consider a Lie group G with H denoting its Cartan subgroup (or maximal 
torus). The required generalized action is given by [3], 

Io = ^2 n s dt trK^g^g, (3.5) 

where g G G and n s are the highest weights defining a unitary irreducible representation of G, 
Kr s \ are the diagonal generators of the Cartan subalgebra H of G. The summation in (|3.5|) 
extends up to the rank r of the algebra, but some of the n's could vanish. 

Given this structure, we see that field theoretic generalization, which would give rise to 
(|3.3|) . appears as J dtdr Y7 s =i^str g~ 1 g, where now n s and g depend on r, while the 
Kr s \ remain as constant elements of the Cartan subalgebra of the group. Thus we take as the 
Lagrange density for our non-Abelian fluid dynamics the formula [3], [Jj 

r 

C = J2i(s) 2 tv K (s)9' 1 D fl g - f(n {1) ,n {2 ), ■ ■ ■ ,n(r)) + C gauge . (3.6) 

Here are a set of Abelian currents; they may be taken to be in the Eckart form = 

n ( s ) u (l)' wriere u/i ' s are four-velocity vectors, with u^u^^ = 1; C gauge governs the gauge 
dynamics of the fields to which the fluid couples. As far as the variational problem of this 
action is concerned, we regard as given by via = y/j^j^s) ■ The space and time 

components of the currents are given by jj\ = (p( s )> P(s) v (s))i P(s) = n (s) n ('s)- ^ n equation 
rt( s ) are the invariant densities for the diagonal directions of the Lie algebra. 
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Comparison with the usual form of the action for an Abelian fluid shows that what we 
have obtained is the non- Abelian analogue of the irrotational part of the flow. In the Abelian 
case, and without the gauge field coupling, equation 1)3.6)1 entails a single contribution, 8 = 1, 
and 2 tr(Kig~ 1 d^g) = —d^O, with vanishing vorticity. In the non- Abelian case, the vorticity 
is nonvanishing. One can nevertheless generalize (j3.fi)) to include the other components of the 
Clebsch parametrized vector (|1.1U|) which couples to This gives the Lagrange density [I], 





r 



£ = ^2j( s ) { 2tlK (s)9 1 D, M g + a^s)} ~ 7(™(1), "(2), • • • , n (r)) + £ gauge, (3.7) 
s=l 

where cw s ) is given by 

= "(s)^/?(s). (3.8) 

For the rank - one group SU (2), with its single Cartan element, the s-sum in ()3.6|) is exhausted 
by a single element. We see that the SU(2) fluid has one component. For higher rank groups, 
the non- Abelian fluid can have up to r components, but fewer — indeed even just a single flow — 
are possible when some of the densities nt s \ vanish. Mathematically, single-component fluids 
are the simplest, but physically it is unclear what kinematical regimes of a quark-gluon plasma, 
for example, would admit or even require such a reduction in flows. 
The covariant derivative in ()3.6)) . namely, 

D„g = d^g + A„g, (3.9) 

involves a dynamical non-Abelian gauge potential = A a JT a whose dynamics is provided 
by C gauge- The first term in C contains the canonical 1-form for the fluid variables in our 
theory and determines their symplectic structure and canonical brackets. We have added 
the Hamiltonian density part, the function /(nm, n (2)i ■ ■ ■ > n (r))i which describes the fluid 
dynamics. The theory is invariant under gauge transformations with group element U , 



9 - U- X g 
A» - U^iA^ + d^U . 



-It A l £t \ TT (3-10) 



where U £ G. One can show that the canonical structure of our theory leads to the charge- 
density algebra ((3~3|> |H, [7]. 

The current that couples to the potential A^ in Q3.7JI reads (matrix rotation) 



r 

JM = 2 Q ^ j (s) ' with = 9K {s) g- 1 . (3.11) 

8=1 

Invariance against arbitrary variation of g ensures that ()3.11|) is covariantly conserved as in 
1)1.11(1 . D^J^ = 0, |5, HI but we also need the conservation of individual [3], 0. This is 
achieved by considering special variations of g of the form 5^g = gKr s \\r s \. These variations 
of g lead to 

d ^{s) = Pis) + V • (v (a) p w ) = 0. (3.12) 

The fluid Wong equation which follows from the conservation of the Abelian and non-Abelian 
currents now reads 

r 

5>£)^V2« = 0. (3.13) 

8=1 
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Varying the individual in (|3.7|) (recall that = yj^3n(s)) produces the non-Abelian 
Bernoulli equations 

2tr Q {s) {D^g)g- 1 =« M /W, (3.14) 

where 

d 

f(S) = 7h7~ f(n { i),n {2) ,...,n (r) ). (3.15) 

un (s) 

The curl of equation ()3.14j) can be cast in the form 

dfl " ° U = 2tr (iD"Q(.))(D v g)9- 1 + Q (S )F^ • (3.16) 

When contracted with j?-. = n^u^y this leaves 

Wfofy (< S )/ (S) ) - "W^/ (S) = -V( S ) 2tr ((D' t Q {a) )(D v g)g- 1 + Q W F^ . (3.17) 

In the single channel case, the right side simplifies, since the sums in (|3.11|) and ()3.13j) contain 
only a single term. Suppressing the subscript (s), we are left from Q3.17|) with, 

nu» d M {u p f{n)) - n d v f (n) = 2 tr J M F» v . (3.18) 

The left side is of the form of the usual Abelian Euler equation [compare (|1.8[) ]: the right side 
describes the non-Abelian Lorentz force acting on the charged fluid. 

In the multichannel case, (|3.11|) an d (|3.13|) are not available because no summation is 
performed on (|3.17jl . and we must remain with that equation - the non-Abelian Euler equation 
for multi-channel flow. 

The matter part of the energy-momentum tensor is given as 



= -r E n Mf {s) - f + E u U u U n (°)f (s) ■ ( 3 - 19 ) 

\s=l / 8=1 

Its divergence 

r 

d ^ v = E { d M'K)) u Uf (s) + n (s) ( u U d » K«)/ (s) ) - dV f {s) ) } ( 3 - 20a ) 

s=l 

is evaluated with (|3.12|) and (|3.13() as 

E^W 2tr ({D^Q {s) ){D^g)g- 1 + Q (s) FA. (3.20b) 

s=l ^ ' 

Since now we are summing over all channels, it follows from (|3.11|) and ()3. 17|> that, 

d^ u = 2tvJ^ u , (3.20c) 
which is canceled by the divergence of the gauge- field energy-momentum tensor. 

d^ auge = -2trJ,F^ (3.21) 
Thus we have conservation of the total energy-momentum tensor. 
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Some simplifications, which lead to a more transparent physical picture, occur if the dy- 
namical potential separates. 



/(ra (1) ,...,ra (r) ) = ^ /(.)("(,)) 



s=l 

f (s) = f{ s) (3-22) 



Then the left side of 1)3.17)1 refers only to variables labeled s, while the right side may be 
rewritten with the help of the non-Abelian Wong equation ()3.13)1 to give 

n M«f.) °i> { u Uf(s))- n (s)9 U f[ s) = 2tr Ij^-j^Uis') [Q{s')F^ + (D^Q {s , ) )(D"g)g- 1 ) 

^ s'/s 

(3.23) 

Thus in the addition to the Lorentz force, there are forces arising from the other channels 
s' ^ s. 

Note that for separated dynamics ([3.22ft . the energy- momentum tensor 1)3.19)1 also sepa- 
rates, 

^ = E % = E {-ST (n {s) f' (s) - / w (n w )) + ^4,} • (3.24) 

8=1 S=l 

but the divergence of each individual 9^ does not vanish. This is just as expected since 
energy can now be exchanged between the different channels and with the gauge field; this 
is also evident from the equation of motion 1)3,23)1 , It is clear that this fluid moves with r 
different velocities vv s ). 

The single-channel Euler equation 1)3.18)1 is expressed in terms of physically relevant quanti- 
ties (currents, chromomagnetic fields); the many-channel equation 1)3.17)1 involves, additionally, 
the gauge group element g. One may simplify that equation by going to special gauge, for 
example g = I, so that the right side of 1)3. 17|) reduces to 

i£ } 2tr ({D^Q is) ){D v g)g- 1 + Q (a) F^\ = f {s) 2tx K {s) {d^A u - fyAJ, (3.25) 

while the Wong equation 1)3.13)) becomes 

r 

£jf s) LV# (s) ]=0. (3.26) 

8=1 

It is interesting that in this gauge the nonlinear terms in disappear. 

Observe that the inclusion of the a, /3-components of the Clebsch parameterizationcurrent, 
or the use of the Lagrangian 1)3.7)1 , does not change the form of the equations of motion for fluids 
when expressed in terms of the velocities and densities. The expressions for these quantities 
in terms of the group parameters and in terms of a, (3 will, of course, be altered. 

We record the nonrelativistic limit of the single-channel Euler equation 1)3.18)1 . (For this 
calculation we restore the velocity of light c, which had been scaled to unity.) For small 
velocities, we may write 



1 

2~ 2 

c,v). (3.27) 
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Further, we take f(n) to be of the form 



/ = nC 2 + V{n) . (3.28) 

With these simplifications, we find that the nonrelativistic limit for the spatial component of 
(|3.18j) gives the Euler equation with a non-Abelian Lorentz force, 

v + v Vv = force + Q a E a + ^ x Q a B a , (3.29) 

where force is the pressure force coming from the potential V (and is therefore Abelian in 
nature), while non-Abelian force terms involve the chromoelectric and chromomagnetic fields. 

£1 = ^, Bi = ~^ k Ff k (3.30) 

It is seen that the non-Abelian fluid moves effectively in a single direction specified by j = pv. 
Nevertheless, it experiences a non-Abelian Lorentz force. 



4 Field-Based Fluid Mechanics 

The theory that we have developed arose from an underlying particle picture of matter. Even 
before committing to definite dynamics (single-channel or multi-channel), very general con- 
siderations on the particle substratum led to the non-Abelian Eckart parameterization (|2.18|) 
and the fluid Wong equation (|2.17j) . For contrast, I now present a derivation of different non- 
Abelian fluid equations, which arise from a field (as opposed to particle) description of the 
fundamental substratum. This is the Madelung representation for the Schrodinger theory: in 
the Abelian case the resulting fluid equations are of the usual Eulerian form. (This is in keep- 
ing with the fact that the Abelian Eckart form is identically true; it requires no committment 
to a specific physical picture.) In the non-Abelian case the equations differ from what has 
been previously derived within a particle description of the underlying matter. 

A field-based realization of the Euler equations for an Abelian fluid is provided by the 
Madelung "hydro dynamical" rewriting of the Schrodinger equation jllj . 

h 2 

ihip = VV (4.1) 

(We consider only the free equation.) Upon presenting the wave function as 

^(t,r) = v 7 ^) e me(t ' r)M , (4.2) 

we find that the imaginary part of (|4.1|) results in the continuity equation for = (p, j), where 
the spatial current j is also the quantum current (h/m) Im^*V"0. When j is written as vp, 
v is identified as V#; the velocity is irrotational. The real part of (|4.1|) gives the Bernoulli 
equation, with a quantum force derived from 

2m 2 ^ 8m 2 p ^ 

The Euler equation follows by taking the gradient of the Bernoulli equation. In this way, we 
arrive again at the conventional irrotational fluid. 
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The story changes if we start from a non-Abelian Schrodinger equation. Again we consider 
the free case with nonrelativistic kinematics. Thus the equation involves a multi-component 
wave function ^f, with 

h 2 

im = V 2 f. (4.4) 

2m 

The color degrees of freedom lead to the conserved non-Abelian current. 

p a = i^T a ^, J a = — Re^ f T a V^ (4.5) 

m 

The singlet current = (p,j) is also conserved. For definiteness and simplicity, we shall 
henceforth assume that the group is SU(2) and that the representation is the fundamental 
one: T a = a a /2i, {T a ,T b } = —5 a b/2. We shall also set the mass m and Planck's constant h 
to unity. The non-Abelian analogue of the Madelung decomposition (|4.2|) is [I], [7], ^21 

V = y/pgu, (4.6) 

where p is the scalar Wty, g is a group element, and u is a constant vector that points in a 
fixed direction [e.g., for SU{2) the two-component spinor u could be taken as u\ = 1, ui = 0, 
then iu^T a u = S^/2]. The singlet density is p, while the singlet current j = Im^*^^ is 

j = vp, v = -iu 1 g~ x Vgu . (4.7) 

With the decomposition (|4.6|) . the color density l|4.7|) becomes 

Pa = QaP, Q a = iu ] g- l T a gu = iR ah u ] T b u= R ab t b /2, (4.8) 

where R a b is in the adjoint representation of the group and the unit vector t a is defined as 
t a /2 = iu^T a u. On the other hand, the color current reads 

Ja = \p Rab u f (T b g~ l Vg + g-'Vg T b ) u , (4.9) 



which with the introduction of 



may be presented as 



g~ x Vg = -2v a T a , (4.10) 

v = v a t a , (4.11) 

Ja = P ~Rab V & • (4.12) 

Unlike the Abelian model, the vorticity is nonvanishing. 

V x v a = e abc v b x v c (4.13) 

A difference between the Madelung approach and the previous particle based one is that 

the color current is not proportional to the singlet current, as in the non-Abelian Eckart 
parameterization. Equation Q4.12JI may be decomposed as 

Ja = Qa/OV + ^RabA (4.14) 
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where the "orthogonal" velocity v" is defined as 



Equation (|4.14f) shows that color current possesses components that are orthogonal to the 
singlet current. 

The decomposition of the SU(2) Schrodinger equation with the parametrization Q4.6JI results 
in two equations: one regains the conservation of the Abelian current and the other is the 
"Bernoulli" equation. 



t a + -V • (pe abc v b t c ) (4.16) 



It is further verified that the covariant conservation of the color current is a consequence of the 
Abelian continuity equation and of (|4.16j) . However, there is no Wong equation because the 
color current is not proportional to the conserved singlet current. Finally, using the identity, 
which follows from the definition (|4.1()jl . 

v a = -\ V [g-HT + e abc v b (g- l g) c , (4.17) 

one can deduce an Euler equation for v a from (|4.16|) . No equation emerges that is analogous 
to the non- Abelian fluid equation (|2.17|) . 

It is clear that within the present approach the fluid color flows in every direction in group 
space, but the mass density is carried by the unique velocity v. This is in contrast to our 
previous approach where all motion is in a single direction or at most in the directions of 
the Cartan elements of the Lie algebra. The difference between the two approaches is best 
seen from a comparison of Lagrangians. For the color Schrodinger theory in the Madelung 
representation 

CSchrodinger = ~ (* f * ~ & *) ~ ~ V** • V* (4.18) 



= ipvt g~ x gu- ipv a • v a - (4.19) 
With u®u^ = 1/2 — 2iK, the free part of the above reads 

^Schrodinger = P 2 tr Kg~ l g - ^p V a ■ v" . (4.20) 

On the other hand, the free part of the Lagrange density (j3,6|) in the nonrelativistic limit, with 
f(n) given by nc 2 = \J p 2 (c 2 — v 2 ) (we restore the velocity of light c), reads 

£° = p 2tr (Kg- l g + Kw ■ g^Vg) - V P 2 (c 2 - v 2 ) , 

« p2tr (Kg- l g + Kvg- l Vg\ - pc 2 + ^pv\ (4.21) 

= p 2 tr Kg~ x g - ^pv 2 - pc 2 , 

where we have used v = — 2tr Kg~ 1 Vg, which follows upon the variation of v, in the next-to- 
last equality above. Thus the canonical 1-form is the same for both models while the difference 
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resides in the velocity dependence of their respective Hamiltonians. Only the singlet v enters 
(|4.21|) while the Madelung construction uses the group vector v a . 

Finally, note that while the Euler equation, which emerges when (|4.16|) and (|4,17|) are 
combined, intricately couples all directions of the fluid velocity v a , it does admit the simple 
solution v a = vt a , with v obeying the Abelian equations that arise from (|4.1|) - (|4.2fl . 
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